VERIFYING ASTEROSEISMICALLY DETERMINED 
PARAMETERS OF KEPLER STARS USING HIPPARCOS 
PARALLAXES: SELF-CONSISTENT STELLAR PROPERTIES 
^ . AND DISTANCES 

V. Silva Aguirre 1 ' 2,3 , L. Casagrande 4 , S. Basu 5 , T. L. Campante 6,7 , W. J. Chaplin 7,1,3 , 
D. Huber 8 , A. Miglio 7 , A. M. Serenelli 9,3 , J. Ballot 10,11 , T. R. Bedding 12,1 , 
J. Christensen-Dalsgaard 1,3 , O. L. Creevey 13 , Y. Elsworth 7 , R. A. Garcia 14,3 , 
R. L. Gilliland 15 , S. Hekker 16,7 , H. Kjeldsen 1 , S. Mathur 17,3 , T. S. Metcalfe 18 , 



< 

M. J. P. F. G. Monteiro 6 , B. Mosser 19 , M. H. Pinsonneault 20,3 , D. Stello 12 , A. Weiss 



Or 

6 



> 
ov 

(N 

00 
O 
(N 



2 

C11U , -TV. VVC1SS 

P. Tenenbaum 21 , J. D. Twicken 21 , and K. Uddin 22 



- 2 - 



ABSTRACT 



1 Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 
120, DK-8000 Aarhus C, Denmark 

2 Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, D-85748 Garching bei Munchen, 
Germany 

3 Kavli Institute for Theoretical Physics, Santa Barbara, CA 93106, USA 

4 Research School of Astronomy and Astrophysics, Mount Stromlo Observatory, The Australian National 
University, ACT 2611, Australia 

department of Astronomy, Yale University, P.O. Box 208101, New Haven, CT 06520-8101, USA 

6 Centro de Astroffsica and Faculdade dc Ciencias, Universidade do Porto, Rua das Estrelas, 4150-762 
Porto, Portugal 

7 School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK 
8 NASA Ames Research Center, Moffett Field, CA 94035, USA 

9 Instituto de Ciencias del Espacio (CSIC-IEEC), Facultad de Ciencies, Campus UAB, E-08193 Bellaterra, 
Spain 

10 CNRS, Institut de Recherche en Astrophysique et Planctologie, 14 avenue Edouard Belin, F-31400 
Toulouse, France 

"Universite de Toulouse, UPS-OMP, IRAP, F-31400 Toulouse, France 

12 Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2006, Australia 

13 Laboratoire Lagrange, UMR 7293, Universite de Nice Sophia-Antipolis, CNRS, Observatoire de la Cote 
dAzur, BP 4229, 06304 Nice Cedcx 4, France 

14 Laboratoire AIM, CEA/DSM-CNRS-Universite Paris Diderot, IRFU/SAp, Centre de Saclay, F-91191 
Gif-sur-Yvette Cedex, France 

15 Center for Exoplanets and Habitable Worlds, The Pennsylvania State University, University Park, PA, 
USA 

16 Astronomical Institute "Anton Pannekoek," University of Amsterdam, Science Park 904, 1098 XH Am- 
sterdam, The Netherlands 

17 High Altitude Observatory, NCAR, P.O. Box 3000, Boulder, CO 80307, USA 

18 Spacc Science Institute, Boulder, CO 80301, USA 

19 LESIA, CNRS, Universite Pierre et Marie Curie, Universite Denis Diderot, Observatoire de Paris, F- 
92195 Meudon Cedex, France 

20 Department of Astronomy, The Ohio State University, Columbus, OH 43210, USA 

21 SETI Institute/NASA Ames Research Center, Moffett Field, CA 94035, USA 

22 0rbital Sciences Corporation/NASA Ames Research Center, Moffett Field, CA 94035, USA 



-3 - 



Accurately determining the properties of stars is of prime importance for 
characterizing stellar populations in our Galaxy. The field of asteroseismology 
has been thought to be particularly successful in such an endeavor for stars in 
different evolutionary stages. However, to fully exploit its potential, robust meth- 
ods for estimating stellar parameters are required and independent verification of 
the results is mandatory. With this purpose, we present a new technique to ob- 
tain stellar properties by coupling asteroseismic analysis with the InfraRed Flux 
Method. By using two global seismic observables and multi-band photometry, 
the technique allows us to obtain masses, radii, effective temperatures, bolomet- 
ric fluxes, and hence distances for field stars in a self-consistent manner. We 
apply our method to 22 solar-like oscillators in the Kepler short-cadence sample, 
that have accurate Hipparcos parallaxes. Our distance determinations agree to 
better than 5%, while measurements of spectroscopic effective temperatures and 
interferometric radii also validate our results. We briefly discuss the potential 
of our technique for stellar population analysis and models of Galactic Chemical 
Evolution. 

Subject headings: asteroseismology — parallaxes — stars: distances — stars: 
fundamental parameters — stars: oscillations 



1. Introduction 



Studying the structure and evolution of the Milky Way requires detailed knowledge of 
the properties of the stellar populations comprising it. In this respect, asteroseismology is a 



powerful too 



Mosser et al 



2006, 



to determine masses and radii of single stars to a high level of p recision (e.g., 



2010 



Michel et al 



Kallinger et all l2010l Metcalfe et all 12 010) 



The CoRoT ( 


Baelin et al. 


2010 




Borucki et al. 


2010. 



Koch et al.ll2010l ) have provided data on stellar oscillations of exquisite quality for thousands 
of stars, encouraging us to carry out a complete stellar census of the observed populations. 

From the thousands of light curves obtained by the space missions, two asteroseismic pa- 



rame ters can be readily extracted (e.g-. lHekker et al.ll2011al . lHuber et al.ll201ll . IChaplin et al. 



20111 ). First, the power spectrum of solar-like oscillators is modulated in frequency by a 
Gaussian-like envelope, where the frequency of maximum power z/ max scales approximately 
with the surface gravity and effective temperature. Second, the near-regular pattern of high 
overtones presents a dominant frequency spacing called the large frequency separation, Au, 
which scales approximately with the square root of the mean stellar density. Applying scal- 
ing relations from solar values, these two asteroseismic observables may be used to estimate 
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stellar properties of large numbers of solar-like oscil l ators, where individual frequencies are 
not available for all targets fe.g.. ISte 



Huber et al.ll201ll ISilva Aguirre et al 



20114 



lo et al. 2008. Basu et al. 2010. Hekker et al. 2011a 



To gain insight about the formation history and evolution of our Galaxy, characteristics 
of stellar populations distributed across it must be accurately known. The best studied 
sample of stars in the Milky Way is the solar neighborhood, where observations and analysis 
over many years have determined some of its key properties, suc h as orbits, kinematics 
and metallicities ( e.g. 



Edvardsson et al. 



1993, 



Reddv et al. 



van Leeuwen 2007 , Feltzing fc Bensby 20081 . Casagrande et al. 2011 ). These data comprise 



2003, Nordstrom et al. 2004 



the basic set of constraints for any model of chemical evolution of the Galaxy. 

Models of Galactic Chemical Evolution are constructed under certain assumptions re- 
garding the physical processes involved in the evolution of our Galaxy, and then calibrated 
against available observations. Those that reproduce them successfully are also used to pre- 
dict other properties of the Galaxy, suc h as abundance gradients across the disk, gas infall 



episodes, and star formation rat es (e.g.. lTinsleylll980l . IChiappini et al.lll997l . iPortinari et al. 



19981 . ISchonrich fc Binneyl 120091 ) . Thus, their predictive power for our galactic history and 
morphology critically depends on how well they can reproduce these observations, most of 
which come from the solar neighborhood sample. Of particular importance among these 
restrictions is the age-metallicity re lation, constructed using stellar isochrones and determi- 



nations of element abundances (e.g-. lEdvardsson et al.lll993l . INordstrom et al.ll2004j ). The ex- 



istence of an age-metallicity relation in the solar neighborh ood is still a subject of debate (see 
Feltzing et al.ll200ll . INordstrom et al.ll2004L lFreemanll2012l . and references therein), and accu- 
rate age determinations are of prime importance to shed some new light in this issue. How- 
ever, the solar neighbo rhood sample used to c onstrain these models is only complete to dis- 
tances of ~50 pc (e.g.. INordstrom et al.ll2004j ). and accurate pr operties of stars further than 
~100 pc are difficult to measure, yet of crucial importance (e.g-. lFreeman fc Bland-Hawthorn 
20021 . ISteinmetz et al.ll2006l llvezic et al.ll2008l ). To extend the sample used as a testbed for 
comparison, we need stellar parameters measured with high accuracy in different regions of 
the Galaxy. 

Asteroseismology can help bridge this gap by providing accurate stellar properties, in- 
cluding distances, for field stars out to several hundred parsecs. These parameters, together 
with effective temperatures, metallicities, and kinematics, should make possible to study 



Midio 


2012a. 


Cheng et al. 


2012, 


Miglio et al. 


2012c) 



Moreover, robust age determinations obtained combining this information with evolutionary 
models will allow construction of the age-metallicity relation of the stellar populations ob- 
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served by CoRoT and Kepler, and so provide tests outside the solar neighborhood of Galactic 
Chemical Evolution models. 

The first comparison between asteroseismically det ermined parameters and predictions 
from Galactic Che mical Evolution mod els was made by iMiglio et al.l (120091 ) for a sample of 
CoRoT red giants. IChaplin et al.l ( 120111 ) used the same technique to obtain masses and radii 
of Kepler main-sequence and subgiant stars, and found a slight but statistically significant 
differ ence between the ob served and synthetic mass distributions. Continuing this line of 
work, IMiglio et al.l ( 12012bl ) obtained distances to CoR oT and Kepler r ed gia nts using bolo- 
metric corrections retrieved from the literature, while ICreevey et al.l ( 120121 ) took a similar 
approach for five Kepler subgiant stars. 

Considering the enormous potential of the results, it is important to verify the tech- 
niques applied in asteroseismic analysis. So far, empirical tests of t he scaling relations 
have used heterogeneous samples and relied on evo lutionary models (see lStello et al.ll2009al . 
Migliol l2012al iBeddingl 1201 ll . iMorel fc Mi grid 120121 and references therein) . In this paper, 
we present a new method to derive stellar parameters in a self-consistent manner, combining 
seismic determinations with the InfraRed Flux Method (IRFM). We compare our results 
with Hipparcos parallaxes, high-resolution spectroscopic temperature determinations, and 
interferometric measurements of angular diameters. We briefly discuss the implications for 
models of Galactic Chemical Evolution and for age determinations of main-sequence stars. 



2. Sample selection and data extraction 



From the more than 500 main-s equence and subgian t stars in which Kepler detected os- 
cillations in its short- cadence mode (|Chaplin et al.ll201ll ). we selected the 32 stars that have 
Hipparcos parallaxes ( Ivan Leeuwenl 120071 ) . We discarded targets with parallax uncertainties 
larger than 20% as we ll as multiple syste ms and incorrect identifications from the Kepler 
Input Catalogue (KIC; iBrown et al.l 120 111 ). This left us with a final sampl e of 22 stars. We 
also retrieved multi-band BtVt and JHK $ photometry from th e Tycho2 ( IH0g et al.ll2000l ) 
and Two Micron All Sky Survey (2MASS; ISkrutskie et al.l 120061 1 catalogs, respectively. We 
have not used the uqriz photo metry from the KIC since it suffers from zero-point uncertain- 



ties ( jPinsonneault et al. 



20121). while Tych o2 data have been extensively tested and tied to 



stellar parameters ( jCasagrande et al.l 120 101 ) . 



Using Kepler data, corrected as described by iGarcia et al.l (120111 ). the global oscillation 
observables z/ max and Av were obtained in two ways. Set A was generate d for all 22 stars 
from the power spectra using the pipeline described by lHuber et al.l (120091 ). while set B was 
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derived from the individu al mode frequencies for the 19 stars for which these are published 
( lAppourchaux et al.ll2012l ). Set B was based on post-survey data only, meaning that at least 
three-month-long time series have been used in the analysis. The analysis of the three stars 
belonging exclusively to set A (namely, KIC 5774694, KIC 5939450 and KIC 10513837) was 
based on survey data, i.e., of one- month- long duration. Set A includes all of set B. 

The way in which seismic parameters were extracted from the frequency lists (set B) 
deserves a brief explanation: i/ max was first estimated by fitting a Gaussian function to the 
envelope of radial (I = 0) mode amplitudes as a function of frequency. This least-squares fit 
was weighted by the uncertainties on the observed amplitudes. Using this value of u maiL1 we 



computed a proxy of Au, Au proxy , using the scaling relation given by IStello et al.l (j2009al ). 
This proxy was then used to construct a Gaussian centered at z/ max with an FWHM of 
AAu proxy , which served as the statistical weight when performing a l east-squares fit to the 
radial frequencies as a function of radial order n to determine Au (see lWhite et al.ll2011al Jbl). 
The choice of such a wide weighting envelope around i/ M v ensures that any f r equency depen 



dence of Au due to acoustic glitches is averaged out (see iMosser et al.ll201ll . iKallinger et al. 



2012|). 



The second step in the selection of a final set of asteroseismic input parameters consisted 
of performing an internal check on the quoted (uncalibrated) uncertainties. The adopted pro- 
cedure makes use of the results obtained for the stars common to both sets. We started by 
computing the rms of the relative residuals in z/ max and Au, i.e., (set B— set A)/set A and 
(set B— set A)/setB, to be subsequently used in the calibration of the quoted uncertainties 
given in sets A and B, respectively. This value of the rms, which can be regarded as an addi- 
tional fractional uncertainty, was then added in quadrature to the uncalibrated uncertainties 
in order to obtain the final (calibrated) uncertainties. 

Since our goal is to perform a uniform analysis, we must use a set of seismic parameters 
obtained from a single method of extraction. Set A has been selected as the final set (after 
error calibration) because it covers all the stars. Fractional differences between results in sets 
A and B are less than 3% for u max and below 0.5% for Au, which are within the uncertainties. 
The procedure used to obtain set A has been extensively tested for consistency with other 
methods of extracting seismic parameters (see lVerner et al.ll201ll and references therein). 



3. Determining stellar parameters 



Our technique for estimating stellar parameters relies on an iterative approach that 
couples asteroseismic analysis with the results of the IRFM. We do this in such a way that 
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the final values of T e g, bolo metric flux, mass and radius are dependent on one another and 
internally consistent. 



3.1. The Direct and Grid-based Methods 



To very good approximation, Au scales as the square root of the mean density (e.g., 
Ulrica Il9861h whil e u max is related to the acoustic cutoff frequency of the atmosphere (e.g., 
Brown et al.l Il99ll . iKjeldsen fc Bedding! Il995l . iBelkacem et al.l l201ll ) . T hese two quantities 



follow scaling relat ions from the accurately known solar parameters (e.g.. lHekker et al.ll2009 
Stello et al.l l2009al ) . which can be written as 



M 



1A 



max,0 



Au 
Au^ 



Au 



T, 



off 



T, 



Here, T f 



eS,Q 



R 

Rq \ 

5777 K, Au Q = 135.1 ±0.1 /iHz and u maX)& 



T, 



efi 



T, 



eff,0 



3/2 



1/2 



(2) 



3 090 ± 30 AtHz are t he observed 



values in the Sun, derived using the same method as set A (IHuber et al.l 120111 ). Provided 
a value of T e g is available, these scaling relations give a determination of stellar mass and 



radius for each star that is independent of evol utionary models (see, e.g.. iMiglio et al.ll2009 



Hekker et al.ll2011bl . ISilva Aguirre et al.ll2011al ). in what has come to be known as the direct 
method. 

Another approach is to include models of stellar evolution when estimating the masses 
and radii. This so-called grid-based method uses evolutionary tracks constructed with a range 
of metallicities a nd searches for a best-fitting model, using Au, Vr^*, T^, and [Fe/Hj as in put 
parameters (e.g., 



Stello et al. 2009b. Basu et al. 2010, Gai et al. 



2011 



Basu et al.1 12012 ). 



Our reference grid of stellar models was c omputed with the Garching Stellar Evolu- 
tion Code (GARST EC, IWeiss fc Schlattll 120081). We have us ed Irwin's equation of state 
( iCassisi et al.l 120031 ). OPAL opacities (llglesias fe Rogers! 119961 ) that were comp lemented at 
l ow te mperatures by those of lFerguson et al. fi2005h . and nuclear reaction rates by lAdelberger et a 
( 1201 ll ). Mode ls for masses below 1.4 M & included microscopic diffusion of helium and met- 
als, following iThoul et al.l (119941 ) . This effect was not considered for masses above 1.4 M & 
since its overall evolutionary impact is small and the code does not include other processes 



that might become relevant in these cases, such as radiative levitation (e.g. iTurcotte et al. 



19981 ). Core and envelope convecti ye overshooting has been included in all models, using 
the exponential decay description of iFreytag et al.1 (119961 ) with a scale factor / = 0.02 and a 
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geometric restriction for small convective cores (IMagic et al.ll2010l ). A mixing length param- 
eter of a = 1.811 from a calibrated solar model and an Eddington T — r relation for stellar 
atmospheres have been adopted. 

The grid spans from 0.8 to 2.0 M & in steps of 0.01 M Q . For the metallicity, we de- 
fined [Fe/H] = at the adopted present-day solar photospheric value of Z/X = 0.0230 



(IGrevesse fc Sauvallll998l ). The initial composition of models was, however, defined in terms 
of that of a calibrated solar model with Z = 0.01876 and Y = 0.26896. Because of gravita- 
tional settling, this implies that the change in [Fe/H] of a solar- metallicity track in 4.57 Gyr 
of evolution is actually ~ +0.06 dex. Initial [Fe/H] values in the grid ra nge from -0.54 to 



+0.3 6 dex in steps of 0.1 dex; the adopted A Y/A Z = 1.4 relation (e.g., ICasagrande et al. 
20071 ) allows for the complete determination of initial composition of models. The grid was 
restricted to models with ages above 0.05 Gyr to avoid degeneracy with pre-main-sequence 
models. The largest time step in the main sequence was constrained to 10 Myr, resulting in 
a very dense grid that comprises, in total, about 4.5 x 10 6 models. 

When applying the grid-based approach, we obtained the Au value of each model us- 
ing fr equencies of individual radi al modes calculated with ADIPLS (jChristensen-Dalsgaard 
20081 ). in the manner described by lWhite et al.l (l2011bl ). On the other hand, v max was always 
computed from the acoustic cutoff frequency relation. In order to find the stellar parameters 



best fi tting the input data, we followed a procedure similar to that presented by iBasu et al. 



(120101 ). Briefly, we produced 10,000 Monte Carlo realizations of sets of input parameters us- 
ing random Gaussian noise around the central observed values, and calculated the likelihood 
for models within 3a of all the observables. Considering only those results with a likelihood 
larger than 95% of the maximum likelihood, we formed the probability distribution and 
assigned the uncertainties to be 34% of either side of the median value. 



As it has been extensively discussed in iGai et al.l (120111 ) and IBasu et al.l (120121 ). it is 
important to consider the dispersion in the grid results arising from the use of different 
evolutionary codes and input physics. For instance, the change in T e g during the main- 
sequence phase due to microscopic diffusio n for stars more massive than ~1.3 M can reach 
approximately 150 K (ITurcotte et al.lll998l ). Since this effect is not taken into account in our 
GARSTEC reference grid for masses above ~1.4 M , when applying the grid-ba sed method 



we als o obtained stellar parameters using the Yale-Yonsei isochro nes presented inlBasu et al. 



(2010 ), the evolutionary tracks from Dotter (IDotter et al.l 12008) and Marigo (jGirardi et al. 



20001 . iMarieo et al.ll2008f ). the YREC set of models presented in lGai et all (120111 ). and evolu- 
tionary tracks constru cted with non-sol ar values of the mixing length parameter of convection 
(the "MLT" set from lBasu et al.ll2012l ). These sets of evolutionary tracks were constructed 



using different microphysics (e.g., equation of state, nuclear reactions, and opacities), as well 
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as assumptions in the metallicity scale and treatment of convection, overshooting and grav- 
itational settling. The final results obtained from the grid-based method encompass these 
uncertainties using an error calibration process described in Section 13.31 



3.2. The InfraRed Flux Method 



The IRFM is arguably one of the most direct and least model-dependent techniques to 
determine effective temperatures in stars. It was originally devised to obtain stellar angular 



diam eters with an accuracy of a few percent (IBlackwell &: Shallislll977l. IBlackwell et al. 



1980). Our analysis is based on the IRFM described by ICasagrande et al.l (120061 12010| ). 



1979 



The basic idea is to recover for each star its bolometric J-*b i (Earth) and infrared 
monochromatic flux J r \ m , both measured at the top of Earth's atmosphere. One must then 
compare their ratio to that obtained from the same quantities defined on a surface element 
of the star, i.e. the bolometric flux o"T e 4 ff and the theoretical surface infrared monochromatic 
flux. For stars hotter than ~ 4200 K the latter quantity is relatively easy to determine be- 
cause the near infrared region is largely dominated by the continuum and depends linearly 
on T e g (Rayleigh- Jeans regime), thus minimizing any dependence on model atmospheres. 
The problem is therefore reduced to a proper derivation of stellar fluxes, which can then be 
rearranged to return the effective temperature. Once J-Boi (Earth) and T e g are both known, 
the limb-darkened angular diameter, 9, is trivially obtained. 

In the adopted implementation, the bolometric flux was recovered using multi-band 
photometry ( Tycho2 B^Vt and 2MASS JHKs), and the flux outside of these bands (i.e., the 
bolometric correction) was estimated using a theoretical model flux at a given T e g, [Fe/H] 
and logg. The infrared monochromatic flux was derived from 2MASS JHKs magnitudes 
only. We used an iterative procedure in T e s to cope with the mildly model-dependent nature 
of the bolometric correction and surface infrared monochromatic flux. For each star, we 



used the ICastelli &: Kuruczl (120041 ) grid of model fluxes, starting with an initial estimate of 
its effective temperature and working at a fixed [Fe/H] and log g until convergence in T e g 
within 1 K was reached. 

The uncertainties stemming from the adopted [Fe/H] and logg were taken into account 
in the error estimate, but their importance is secondary at this stage since the IRF M has 
been shown to depend only loosely on those parameters (see ICasagrande et al.l 120061 ) . This 
makes the technique superior to most spectroscopic methods for determining T e g -provided 
reddening is known- since the effects of T e ff, logg and [Fe/H] on the latter are usually strongly 
coupled and the model dependence is much more important. The metallicity adopted for 



each star and the coupling of the IRFM with asteroseismic gravities will be discussed in the 
next section, together with reddening effects. 



3.3. Iterations and Error Determination 

Table 1: Input parameters for the 22-star Sample. See text for details. 



KIC ID 


HIP 


"max (MHz) 


Av (/iHz) 


[Fe/H] 


E(B - V) 


3632418 


94112 


1144 ± 31 


60.8 ± 0.2 


-0.01 


0.024 


3733735 


94071 


2145 ± 61 


92.3 ± 0.3 


-0.10 


0.025 


4914923 


94734 


1887 ± 181 


88.7 ± 0.3 


0.17 


0.018 


5371516 


96528 


1018 ± 33 


55.4 ± 0.2 


0.13 


0.020 


5774694 


93657 


3442 ± 274 


140.2 ± 4.0 


0.01 


0.000 


5939450 


92771 


605 ± 25 


30.5 ± 2.4 


-0.01 


0.020 


6106415 


93427 


2219 ± 60 


104.3 ± 0.3 


-0.06 


0.000 


6225718 


97527 


2338 ± 66 


105.8 ± 0.3 


-0.15 


0.010 


7747078 


94918 


946 ± 26 


54.0 ± 0.2 


-0.26 


0.018 


7940546 


92615 


1081 ± 34 


58.9 ± 0.2 


-0.04 


0.010 


8006161 


91949 


3570 ± 96 


149.3 ± 0.4 


0.34 


0.000 


8228742 


95098 


1175 ± 34 


62.1 ± 0.2 


-0.14 


0.025 


8751420 


95362 


571 ± 15 


34.6 ± 0.1 


-0.20 


0.010 


9139151 


92961 


2695 ± 74 


117.3 ± 0.3 


0.15 


0.012 


9139163 


92962 


1685 ± 45 


81.1 ± 0.2 


0.15 


0.012 


9206432 


93607 


1859 ± 50 


84.7 ± 0.3 


0.23 


0.013 


10068307 


94675 


976 ± 35 


54.0 ± 0.2 


-0.13 


0.014 


10162436 


97992 


1016 ± 28 


55.8 ± 0.2 


-0.08 


0.023 


10454113 


92983 


2310 ± 68 


105.1 ± 0.3 


-0.06 


0.011 


10513837 


91841 


191 ± 7 


14.6 ± 0.2 


0.15 


0.026 


11253226 


97071 


1669 ± 45 


77.0 ± 0.2 


-0.03 


0.011 


12258514 


95568 


1499 ± 40 


75.0 ± 0.2 


0.13 


0.015 



As described in Section 13. 1[ the asteroseismic methods provide a mass and radius based 
on an input T eff value (and [Fe/H] for the grid-based case). On the other hand, the IRFM 
gives T e ff and the bolometric flux at a given input log g and [Fe/H]. In order to determine 
a unique set of stellar parameters for each star, we iterated the two methods in a consistent 
way, using both the dir ect and grid-based approac h. A simplified version of this technique 
was first introduced by ISilva Aguirre et al.l (j2011al ). 
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We started by calculating sets of IRFM effective temperatures for each star at fixed 

log g = 2.0 5.0 in steps of 0.5 dex; this translates into T eff changes of less than 1% for each 

log g step. The metallicity of the targets must be given as an input, and we have considered 
them in the following order of prefer ence, according to availa bility: the latest revision of the 
Gene va-Copenhag e n Sur vey (GCS; ICasagrande et al.l 120111 ). spectroscopic determinations 



from lBruntt et al.l ( l2012f ). or the value given in the KIC increased by 0.18 dex. The latter is 
the offset found between GCS and the KIC for the 1 1 stars common in our sample, and is 
similar to the +0.21 dex offset found by lBruntt et al.l (120121 ). 



Drimmel et a . 


(2003) 


Mielio et al. 


(2 


012bh. 



These were obtained after an iter- 



dependent extinction values from 
ation in distance as described by 
was available we assigned E(B — V) = 0.0. We list in Tabled] the input parameters used in 
our analysis. 

The procedure applying the direct method works as follows. Using log g determinations 
from the KIC as an initial guess, we interpolated in gravity and computed T e g from the 
IRFM results. This T e g value, together with i/ max and Az/, was fed to the scaling relations to 
obtain a mass, radius, and thus logg. Interpolating again in gravity gave an updated value 
of T e ff, and the procedure was repeated until convergence in logg and T e s was reached. 

We obtained la uncertainties of the parameters during the iterations. Uncertainties in 
the seismic observables were taken into account, as well as variations in the T e g determi- 
nations arising from different photometric filters and logg determinations. The results are 
affected by the assumed value of extinction, and are mildly dependent on the metallicity 
considered. To account for possible errors in reddening and composition, we have also com- 
puted sets of results at log g = 3.5, one increasing E(B — V) by +0.01 (the decreasing case 
is essentially symmetric), and another one changing the metallicity by ±0.1 dex. Moreover, 
a Monte-Carlo simulation was run to estimate the uncertainties in T e s from random pho- 
tometric errors. Finally, we added an extra 20 K to the error budget to account for the 
uncertainty in the zero-point of the temperature scale. 

The analysis was repeated using the grids mentioned in Section 13.11 to determine mass, 
radius, and logg values at each iteration. In all cases we used as input values the seismic 
observables and metallicities described above, considering an uncertainty in composition of 
±0.1 dex consistent with what was applied for the IRFM. The final set of stellar parameters 
from the grid-based method and their corresponding uncertainties were obtained in the same 
manner as the seismic input parameters (cf. Section [2]): we adopted the GARSTEC grid as 
the reference and performed an error calibration by computing the rms of the relative resid- 
uals in mass, radius, and gravity, and adding them in quadrature to the original GARSTEC 
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uncertainties. 



4. Results 

The procedure outlined in Section l3~3l provided final values of T e g, mass, radius, log g, 
and, as mentioned in Section 1X2"! the associated bolometric flux and 9. Their corresponding 
la uncertainties were also obtained during the iterations. It is important to mention that 
the log g values determined via the direct and grid-based method agree better than 0.03 dex, 
implying that their T c g- values are also in agreement within 3 K. Using the asteroseismic 
radius and 6, it is straightforward to estimate the distance: 

dseis = C —j— , (3) 

where C is the conversion factor to parsecs. In this manner, we determined asteroseismic 
distances for our 22 sample targets. 

In Figure [1] we compare our distances with those obtained from Hipparcos parallax 
measurements. Note that, as described in Section 13.11 seismic radii determinations can be 
obtained by either the direct or grid-based method. The agreement is excellent, particularly 
for the close-by targets, boosting our confidence on the asteroseismic parameters and the 
robustness of our technique. 

There is one target that clearly deviates from the one-to-one relation in the results 
obtained via the direct method, shown with a red diamond in Figure [TJ Not surprisingly, 
this star has the largest fractional error in Av (~8%, compared to the less than ~2% of the 
rest of the targets, see Equation [2]). The large uncertainty in this parameter dominates the 
radius error budget, and thus propagates to the estimated distance. For this particular case, 
the grid-method determined the radius of the star with a much higher accuracy, a result 
that is confirmed by the better agreement of its distance to that obtained from parallaxes. 
The fact that the grid-based approach restricts combinations of Av, z/ max and [Fe/H] to a 
narrow range of possible T e g values seems to be sufficient to deal with large uncert ainties in 



one measurement, provided the others are accurately determined (IGai et al.ll201lf ) 



For stars with accurate seismic measurements, the direct method provides results as 
reliable as the grid-based one. The weighted mean difference (Hipp. — Seis.) is 2.1% ± 1.8% 
for the direct method, while for the grid-based case is 2.4% ± 1.5%. Removing the outlier 
from the sample changes the average differences to 2.3%±1.8% and 3.1%±1.6%, respectively. 
One important factor to take into consideration is extinction. From Figure [TJ we see that 
the uncertainties in asteroseismic distances seem to increase with distance. This points out 
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Fig. 1. — Upper panels: comparison of Hipparcos distances with those obtained via the 
seismic method. Lower panels: fractional difference (Hipp. — Seis.) between both determi- 
nations. Gray solid line shows the one-to-one correspondence, while the black dashed and 
dotted lines represent the weighted average difference and standard deviation, respectively. 
The red diamond in all panels is the star with the largest fractional error in Av. See the 
text for details. 

to reddening as the cause, since the error in the seismic distance determinations should be 
comparable for star with similar uncertainties in the global seismic parameters. Analysis of 
the error budget shows that reddening becomes the major contributor as distance increases, 
most likely due to the use of distance-dependent integrated maps of extinction. We discuss 
this further in Section [5j 

As described in Section |2j the seismic input parameters were determined using two 
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different methods, and the results shown in Figured] are those from set A (pipeline processing 
of the power spectrum). We have tested the impact on the distance determinations of using 
instead the results from individual frequency lists (set B): for the 19 targets common to 
both sets, containing stars up to ~170 pc away, the weighted mean difference between the 
distances of Set A and Set B is below 0.5%. 



All but one of our targets were included in the spectroscopic analysis made by lBruntt et al 



( 120121 ) . Those authors obtained effective temperatures via the excitation balance of Fe I lines, 
using a fixed log g value in their analysis as determined by asteroseismology. In Figure [2] we 
compare their T e g values with ours and find excellent agreement. Individual fractional differ- 
ences are all below 2%, while the weighted mean difference (Spec. — Seis.) is —0.8% ± 0.4%. 
Th is level of agreemen t is particularly impressive considering that the uncertainties quoted 
by iBruntt et all ( 120121 ) are of 70 K for all the targets. However, there is a possible system- 
atic offset to lower values in the spectroscopic temperatures compared to ours. Although 
the reason for this behavior is far from evident and goes beyond the scope of this paper, 
we mention that the effective temperatures determined in this work are supported by the 
analysis of hydrogen line profiles (M. Bergemann, private communication). The Balmer-line 
T e ff scale is known to b e warmer than that given by the excitation balance of Fe I lines (see 
Casagrande et al.l I2010L Figure 12), making it more consistent with the IRFM (M. Berge- 
mann et al. 2012, in preparation). 

Another verification of our technique comes from comparing our derived angular diam- 
eters with results from interferometry. Four stars in our sample have been o b served with 
the PAVO/CHARA long-baseline interferometer, as described by iHuber et al.l (120121 ) . The 
residual mean value between our a ngular diameters an d the interferometric ones is below 1%, 
consistent with the value found by IHuber et al.l (120121 ) for a larger sample. T hese results are 



also c ompared to the ones derived using the surface brightness technique of iKervella et al. 



( 120041 ). obtaining an equally g ood level of agreem ent and further confirming the robustness 
of our method (see Figure 4 in lHuber et al.ll2012l ). 



The results outlined above clearly show that our method provides accurate T e g, radii, 
angular diameters, and therefore distances for the sample of stars considered. In order to 
assess the level of accuracy of these asteroseismic distances, we separated the sample into 
three bins according to the uncertainty in the parallax measurements. There is a natural 
correlation between distance and quality of the parallaxes, where stars with smaller uncer- 
tainties are usually closer and therefore less affected by interstellar extinction. For those 
stars with parallax determined better than 5%, the rms between our grid-based distances 
and the Hipparcos ones is below 4.7%. Moreover, when considering the bin of uncertainties 
between 5 and 10%, as well as the bin with parallax errors between 10 and 20%, the rms of 
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Fig. 2. — Upper panel: comparison of effecti ve temperatur e s usin g the asteroseismic method 
with those obtained from spectroscopy by iBruntt et al.l (120121 ) . Lower panel: fractional 
difference (Spec. — Seis.) between both determinations. Solid lines shows the one-to-one 
correspondence, while the dashed and dotted lines represent the weighted average difference 
and standard deviation, respectively. 



the relative residuals is also smaller than the typical Hipparcos error. Thus, provided extinc- 
tion is properly taken into account, our distance determinations can be considered accurate 
to 5%. 

We present in Table [2] the parameters obtained with the grid-based method. Note 
that, as mentioned in Section IXTj the grid-based approach returns a distribution probability 
function for each parameter, and so the uncertainties in the mass, radius, and gravity are 
asymmetric. 



- 16 - 



Table 2: Stellar parameters derived using the grid-based method for the 22-star sample. 
Errors on the parameters come from probability distribution functions (see Sects. l3TTl and l3. 31 
for details) 



KIC ID HIP T cff (K) M/M Q R/Rq log g F bol (10~ 8 org s~ 1 cm~ 2 ) 9 (mas) d (pc) 



3632418 


94112 


6286 ± 


70 


1 


396±° 


074 


1 


911+° 


025 


4 


0181° 


003 


1 


402 


± 





074 





104 


± 





006 


108±4 


3733735 


94071 


6824 ± 


131 


1 


454±° 


028 


1 


427^55 


019 


4 


284155 


008 


1 


251 


± 





070 





132 


± 





006 


101±5 


4914923 


94734 


5828 ± 


56 


1 


174+" 

— 


058 
6 2 


1 


408+55 

— 


022 
023 


4 


209+5! 
— u 


003 
004 





456 


± 





023 





109 


± 





003 


120±4 


5371516 


96528 


6360 ± 


115 


1 


468±£ 


063 


2 


066155 


022 


3 


97315! 


007 


1 


167 


± 





062 





146 


± 





007 


131±6 


5774694 


93657 


5911 ± 


106 


1 


079±° 


041 


1 


000+.55 


015 


4 


46715! 


007 


1 


213 


± 





07(1 





173 


± 





008 


54±3 


5939450 


92771 


6380 ± 


70 


1 


625±£ 


134 


2 


916155 


062 


3 


72615! 


008 


3 


207 


± 





169 





241 


± 





008 


113±4 


6106415 


93427 


6061 ± 


89 


1 


110+° 


036 


1 


2401° 


018 


4 


29615! 


004 


3 


491 


± 





024 





279 


± 





008 


41±1 


6225718 


97527 


6338 ± 


88 


1 


209+o 


037 
034 


1 


2561° 


014 
014 


4 


32215! 


004 
004 


2 


664 


± 





291 





223 


± 





014 


52±3 


7747078 


94918 


5856 ± 


112 


1 


135i° 


086 
088 


1 


9521° 


039 
039 


3 


910l° 


004 
004 





452 


± 





09!) 





107 


± 





012 


169±20 


7940546 


92615 


6287 ± 


74 


1 


380+° 


065 
104 


1 


9441° 


024 
031 


3 


99615! 


006 
012 


2 


913 


± 





039 





236 


± 





00(1 


76±2 


8006161 


91949 


5355 ± 


107 





959+. » 


035 
037 





927155 


014 
014 


4 


4841° 


004 
004 


2 


879 


± 





153 





324 


± 





016 


27±1 


8228742 


95098 


6130 ± 


107 


1 


308+° 


062 
060 


1 


8551° 


027 
027 


4 


017l° 


005 
004 





457 


± 





156 





099 


± 





017 


175±31 


8751420 


95362 


5243 ± 


162 


1 


285j;° 


082 
096 


2 


7221° 


048 
057 


3 


6741° 


006 
005 


4 


898 


± 





368 





441 


± 





032 


57±4 


9139151 


92961 


6141 ± 


114 


1 


21S±° 


046 
046 


1 


1781° 


018 
018 


4 


380l55 


004 
004 





527 


± 





017 





106 


± 





004 


104±4 


9139163 


92962 


6525 ± 


111 


1 


405±° 


034 
027 


1 


"it" 


010 
010 


4 


1951° 


004 
004 


1 


225 


± 





031 





142 


± 





005 


103±4 


9206432 


93607 


6614 ± 


135 


1 


4H2±° 


044 
044 


1 


5441° 


015 
015 


4 


23ll° 


005 
005 





605 


± 





0(14 





097 


± 





007 


147±10 


10068307 


94675 


6197 ± 


97 


1 


366+J5 


062 
071 


2 


060+g 


028 
033 


3 


9431° 


003 
004 


1 


401 


± 





018 





169 


± 





005 


114±4 


10162436 


97992 


6245 ± 


110 


1 


3651° 


061 
068 


2 


015l° 


025 
027 


3 


96ll° 


006 
004 





947 


± 





075 





137 


± 





007 


137±7 


10454113 


92983 


6134 ± 


113 


1 


165+° 


045 
045 


1 


2511° 


017 
017 


4 


30915! 


005 
004 





924 


± 





050 





140 


± 





006 


83±4 


10513837 


91841 


4955 ± 


95 


1 


290^55 


072 
076 


4 


788155 


083 
103 


3 


18615! 


006 
007 





585 


± 





043 





171 


± 





009 


261±15 


11253226 


97071 


6715 ± 


97 


1 


4581° 


032 
034 


1 


628155 


017 
018 


4 


1761° 


006 
004 


1 


095 


± 





027 





127 


± 





004 


119±4 


12258514 


95568 


6064 ± 


121 


1 


302+J5 


078 
084 


1 


630155 


029 
031 


4 


1271° 


004 
005 


1 


532 


± 





063 





184 


± 





008 


82±4 



Our parameters can b e compared to othe r stud ies where the same stellar properties 
were determined. R ecently. IPinsonneault et al.l (120121 ) provided color-temperature relations 



consistent with the ICasagrande et al 



scale using the available Sloan Digital Sky 
Survey (SDSS) griz photometry from the KIC. Fourteen stars from our sample are present 
in their catalog, and comparison of the obtained T e g values is shown in Figur^3l The weighted 
mean differe nce (SPSS — Seis.) i s —0.1 % ±0.6%, indicating that the corrected temperatures 
provided by IPinsonneault et al.l (120121 ) are indeed on a scale consistent with ours and can be 



used for studies of field stars when this photometry is available. 



Mathur et al.l (120121 ) performed a detailed modeling of several Kepler targets, including 



six stars from our sample, using individual frequency lists obtained from one-month long 
observations. The input metallicities and effec tive temperatures u sed in that study were in 
some cases different from ours (see Table 3 in iMathur et al.ll2012l ). Comparing the results 
reveals good overall agreement in the obtained radii, with two targets showing what appears 
to be a slight radius underestimation in their determinations with respect to ours. A more 
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Fig. 3. — Upper pane l: comparison of effective t emperatures using the asteroseismic method 
with those derived by iPinsonneault et all ( 120121 ). Lower panel: fractional difference (SDSS — 
Seis.) between both determinations. Solid lines shows the one-to-one correspondence, while 
the dashed and dotted lines represent the weighted average difference and standard deviation, 
respectively. 



thorough investigation of this issue will be made when detailed modeling using longer time 
series is performed (T. M. Metcalfe et al. 2012, in preparation). 

Although our sample of stars only represents a small fraction of the total short-cadence 
Kepler sample, they cover a wide range in metallicity, T e g and logg. In Figure H] we present a 
log g — T c g- diagram, where the 22 targets have been placed using the parameters derived with 
the grid-based method. Also plotted are stellar evolution tracks from the GARSTEC grid, 
at masses and metallicities compatible with those given in Tables [1] and |2j Our targets are 
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Fig. 4. — Position in the log g — T e g- plane of the 22-star sample as derived via grid-based 
approach and the IRFM. Also plotted are stellar evolutionary tracks from the GARSTEC 
grids for masses between 0.9-1.4 M® at two different initial metallicities: [Fe/H] = +0.35 
(dotted blue lines) and [Fe/H] = —0.24 (dashed red lines). 

distributed in different evolutionary stages, from the early main sequence to the beginning 
of the red giant branch. Thus, we have tested the accuracy of our method in stars showing 
a wide variety of masses, ages, compositions, and energy transport mechanisms. 

We find good agreement in the radii and bolometric fluxes (and hence distances) esti- 
mated by the direct and grid-based methods to the parallax data, which again supports our 
previous statement that the method has only a modest sensitivity to the stellar composition 
(see Section I3T2~|) . Most importantly, it shows that we can extend the implementation of our 
method to stars evolved beyond the subgiant branch. This will be addressed in an upcoming 
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publication. 

A particularly interesting case is that of KIC 10513837, the most evolved and most 
distant star in the sample. We estimated its composition to be [Fe/H] = 0.15 by adding 
0.18 dex to the KIC value (see Section [3T31 and Table, [p. However, its position in Figure 0] is 
instead compatible with sub- so lar metallicity. Spectroscopic analysis of this target made by 



Molenda-Zakowicz et al.l (120081 ) found a value of [Fe/H] = —0.07 and an effective temperature 
consistent with our determination, further confirming the mild sensitivity of our method to 
composition. 

In a natural continuation of this work, we applied our procedure using the direct method 
to the complete short-cadence Kepler sample and derived consistent parameters, including 
distances, for all these stars. The 565 targets considered are predominantly main-sequence 
and subgiant stars, with a handful of red giants also present in the sample. All the targets 
have Tycho2 photometry available, and we have used metallicites from the KIC increased 
by 0.18 dex. To account for the uncertainties in composition, we have computed IRFM sets 
of T e ff varying the metallicity by ±0.3 dex (see Section I3~3"|) . The seismic input parameters 
were computed as for set A, described in Section [2J Extinction values were obtained from 



Drimmel et al.l ( 120031 ) after an iteration in distance. These preliminary determinations will 
be compared with results from several pipelines in a forthcoming publication (W. J. Chaplin 
et al. 2012, in preparation). In Figure [5] we show a histogram with the obtained distance dis- 
tribution, where we have also plotted the distribution of our 22 sample stars for comparison. 



The resulting distance distribution shows that we can use Kepler data to probe popu- 
lations of main-sequence and subgiant field stars as far as 1 kpc from the Sun. However, our 
distance determin ations for this sample can be s l ightly undermined by faulty metallicities 



from the KIC (see iMolenda-Zakowicz et al.ll201ll iBruntt et al.ll201ll for a discussion). We 



have estimated [Fe/H] for these stars as consistently as possible with the available data. The 
latest revision of th e GCS is built upon 150 stars with spectroscopic measurements agreeing 
with our T eff scale (jCasagrande et al.ll201ll ). Therefore, on average the +0.18 dex correction 
we find makes the metallicity scale of the KIC consistent with the underlying scale we 



are de riving here. On the other hand, the good agreement with the results of IBruntt et al. 



(120121 ) implies that also the latter metallicity scale is broadly on the same zero point that 
we adopt (and in fact, a similar +0.21 dex correction is found). Thus, although the most 
self-consistent approach would be to iterate our T e g- also for deriving metallicities, on average 
this condition is fulfilled. 



As a final test of the impact on distances of incorrect metallicities, we have done calcu- 
lations considering different assumptions for the composition of the full sample: using the 
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Fig. 5. — Derived distance distributions using IRFM + asteroseismology. Red shaded region 
shows the results for the 22-sample stars, while the unshaded one depicts the results for all 
main-sequence and subgiant stars with available Kepler asteroseismic data. 



[Fe/H] values as provided in the KIC (i.e., not including the +0.18 dex offset), and also 
considering a fixed mean metallicity for all the targets of [Fe/H] = —0.2 dex. The distance 
distribution obtained remained practically unchanged, reinforcing the notion of the mild 
metallicity dependence of our method. 
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5. Conclusions 

Determining accurate stellar parameters is crucially important for detailed studies of 
individual stars, as well as for characterizing stellar populations in the Milky Way. The 
asteroseismic revolution produced by the CoRoT and Kepler missions requires robust tech- 
niques to exploit fully the potential of the data and provide the community with the building 
blocks for ensemble analysis. 

Using oscillation data and multi-band photometry, we have presented a new method to 
derive stellar parameters, combining the IRFM with asteroseismic analysis. The novelty of 
our approach is that it allows us to obtain radius, mass, T e g, and bolometric flux for individual 
targets in a self-consistent manner. This naturally results in direct determinations of angular 
diameters and distances without resorting to parallax information, further enhancing the 
capabilities of our technique. 

Two asteroseismic methods were applied to the available data, one based entirely on 
scaling relations and the other one using grids of pre-calculated models. When accurate 
seismic data are available, comparison of our distance results with those from Hipparcos 
parallaxes shows an overall agreement better than 4%, regardless of the asteroseismic method 
employed. Furthermore, the obtained T e g values show a mean difference below 1% when 
compared to results from high-resolution spectroscopy. We have also compared our calculated 
angular diameters with those measured by long-baseline interferometry and found agreement 
within 5%. This provides verification of our radii, T e g values and bolometric fluxes to an 
excellent level of accuracy. 

Despite the encouraging results, systematics can arise from faulty determinations of 
reddening values and metallicities. In Section [3731 we described how the effects of metallicity 
and reddening, as well as their uncertainties, were taken into account in the calculations. 
For our determinations to be completely self-consistent, we must be able to determine those 
parameters from a single set of data using the IRFM. An observational campaign is currently 
under way to obtain Stromgren photometry of Kepler s tars that will provide a h omogeneous 



set of values for [Fe/H] and extinction, as described by ICasagrande et al.l (120 Iff ). 



For most of the stellar parameters included in our verifications, both asteroseismic 
methods produce equally good results. However, it should be kept in mind that the direct 
method can be significantly biased when large uncertainties in the seismic input parameters 
exist. Moreover, scaling relations are likely to have a different dependance on effective 
temperature bey ond the main-sequen ce phase, as suggested by comparison to evolutionary 



calculations (see IWhite et al.l 1201 lbl ). The restrictions imposed by metallicity and by the 



theory of stellar evolution help to cope better with large errors in seismic data, and the use 
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of the grid-based analysis in these cases is therefore recommended . 

However, to take full advantage of the available parameters, asteroseismology must pro- 
vide masses with a comparable level of accuracy. It is important to note that results on 
masses from the direct method for values above ~ 1.5 M© can deviate significantly from 
those obtained using the grid-based approach. In fact, dif ferences of more than ~30% are 
not unusual in these cases. Using different grids of models, iGai et al.l (120111 ) found that the 
fact that the direct method does not explicitly take metallicity into account could undermine 
its mass determinations. A thorough comparison of different grid-based techniques with the 
direct method is beyond the scope of this paper and will be presented in an upcoming publi- 
cation (W. J. Chaplin et al. 2012, in preparation). Another method to obtain asteroseismic 
mass es is via detailed modeling of targets, aim ing; at fitting the list of individual frequencies 
(e.g. iMetcalfe et al.l |2010L iMathur et al.ll2012[ ). This approach provides mass estimates with 



a high level of precision and, in principle, also with high accuracy. Regardless of the consid- 
ered technique, one must keep in mind that verification of asteroseismic mass determinations 
in general is still needed. 



(Miglio 


2012a. 


Miglio et al. 


2012c) 



Com- 
bining this information with evolutionary models can lead to an age-metallicity relation, 
opening the possibility of testing models of Galactic Chemical Evolution in stars outside th e 



solar neighborhood (e.g., IChiappini et al.lll997l ISchonrich &: Binneyl 120091 iFreemanl 120121 ) . 



Applying our method to the complete short-cadence Kepler sample reveals that we can probe 
stars as far as 1 kpc from our Sun, making this set of main-sequence and subgiant stars ex- 
tremely interesting for population studies. Although much greater distances can be probed 
by analyzing oscillations in giants, the ages of these stars are m ostly determined by their 
main-sequence lifetime (e.g.. ISalaris et al.ll2002l iBasu et al.l 120111 ). Thus, the short-cadence 
sample is of key importance for helping to calibrate mass-age relationships of red giants and 
correctly characterize their populations. 

A substantial number of the Kepler main-sequence and subgiant targets have been 



obser ved long enough to obtain individual frequency determinations (lAppourchaux et al. 



2012h . Detailed modeling of these stars, particularly using frequency combinations(e.g. 
Roxburgh fc Vorontsov[2003. be Meulenaer et al. 2010 ) and modes of mixed character (e.g., 



Deheuvels fc Michel 



201 ll . iBenomar et al.ll2012l ). can put tighter constraints on their masses 



and ages, providing anchor points for ensemble studies. In fact, certain combinations of 
frequencies can be used to probe the remaining central hydrogen conten t in stars (e.g., 
Christensen-Dalsgaardlll988f ) . the existence and size of a convective core (e.g.. lSilva Aguirre et al 
201 lbl . and references therein), and the position of the convective envelope and helium surface 
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abundance (e.g.. IChristensen-Dalsgaard et al.lll99ll lAntia fc Basulll994l ). These techniques 
are currently being applied to several stars in the sample (e.g., S. Deheuvels et al. 2012, 
in preparation, A. Mazumdar et al. 2012, in preparation, V. Silva Aguirre et al. 2012, in 
preparation) and should help us obtain masses with higher accuracy and determine more 
robust differential ages. 
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